Lattice QCD signal for a bottom-bottom tetraquark 
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Utilizing lattice QCD results for the potential of two static antiquarks and two dynamical quarks 
as well as quark model techniques for the dynamics of two heavy antiquarks in a cloud of two light 
quarks, we are provided with an accurate framework for the study of possibly existing heavy-heavy- 
light-light tetraquarks. Among the possible quantum numbers of such a system, we find binding in 
only one channel, the scalar isosinglet. Solving the Schrodinger equation for the displacement of the 
heavy antiquarks and taking systematic errors into account, we find an antibottom-antibottom-light- 
light bound state with a confidence level of around l.Sa . . . 3.0a and binding energy of approximately 
30MeV...57MeV. 

PACS numbers: 12.38. Gc, 13.75.Lb, 14.40.Rt, 14.65.Fy. 



Experimentally exotic hadrons have been searched for 
many years, because as soon as quarks were proposed 
in the sixties, it became clear that systems more com- 
plex than mesons and baryons could possibly exist. How- 
ever, exotic hadrons are very elusive systems. They are 
much harder to observe experimentally, to understand 
by theoretical model calculations, and to simulate by 
means of lattice QCD than the conventional mesons and 
baryons. To confirm the existence or non-existence of ex- 
otic hadrons still remains an important problem in QCD. 

A frequently discussed exotic multiquark is the 
tetraquark. It was already proposed in the seventies [l| 
as a bound state formed by two quarks and two anti- 
quarks. There arc several hadronic resonances, which 
are tetraquark candidates, e.g. a, k, D* q or D sl H. The 
most recent tetraquark candidate has been claimed by 
the Belle Collaboration Q, observing in five different 
T(5S) decay channels two new charged bottomonium res- 
onances Zf, with masses 10610 MeV and 10650 MeV and 
narrow widths of the order of 15 MeV, where the charge 
can only come form the presence of a light quark and a 
light antiquark. However, the tetraquark nature of these 
resonance is disputed [j. 

Notice that the experimental particle physics collab- 
orations arc technically improving rather impressively. 
One decade ago SELEX at FNAL already studied the 
doubly charmed baryons. Presently, BELLE at KEK, 
CDF and D0 at FNAL, and LHCb at CERN have al- 
ready observed bottom hadrons. Thus, they may pos- 
sibly be able to search for charm-bottom and bottom- 
bottom tetraquarks, if a strong evidence is presented by 
theoretical calculations. 

Tetraquark studies face several difficulties. (1) Mesons 
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and baryons only decay strongly, when the confining 
string breaks, a quark antiquark pair is created and and 
cither two mesons or a meson and a baryon are formed. 
In contrast to that, tetraquarks are directly open to 
meson-meson decay. (2) Moreover, tetraquarks are rela- 
tivistic four-body systems, which are highly complex few- 
body systems. (3) And on the top of these technical dif- 
ficulties no model (since the onset of QCD [l[ up to the 
present) seems to be sufficiently well calibrated to ad- 
dress multiquark binding: different quark models, even 
when producing similar meson or baryon spectra, usually 
differ significantly in predictions for tetraquarks. 

A strategy to avoid many technical difficulties asso- 
ciated with tetraquark studies consists in searching for 
bound states rather than for resonances, which is e.g. 
appropriate, when two heavy antiquarks (or equivalently 
two heavy quarks) are involved. This strategy was al- 
ready identified in the eighties On the one hand, it 
is plausible that any light tetraquark can only be treated 
as a resonance, because it couples to two-meson chan- 
nels with identical quantum numbers, where at least one 
of the mesons is a pion or kaon. Due to chiral sym- 
metry, these are very light mesons and, consequently, 
the tetraquark will have a rather light open decay chan- 
nel. On the other hand the presence of two heavy anti- 
quarks is sufficient to force all involved hadrons, i.e. the 
tetraquark and any two-meson channels with the same 
quantum numbers, to escape chiral symmetry. More- 
over, the heavy antiquarks also reduce the technical com- 
plexity of the four quark system, since heavy antiquarks 
are non-relativistic. This in turn allows for the Born- 
Oppenheimer approximation: for the light quarks the 
heavy antiquarks can be regarded as static color charges; 
once the light diquark energy is determined, it can be uti- 
lized as an effective potential for the heavy antiquarks. 

The only remaining theoretical problem is to obtain 
the correct effective potential, which has been studied by 
lattice QCD methods mainly in the quenched approxima- 
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tion (cf. e.g. @, 0] and references therein). Only recently 
computations of this potential with dynamical sea quarks 
have been performed [§-[l(| • Here we use such dynamical 
results obtained with a comprehensive set of four quark 
operators While for the static antiquarks the 

only relevant variable is their separation, the two light 
u/d quarks can be combined in 8 x 8 = 64 different ways 
(each light quark has two isospin, two spin and two parity 
degrees of freedom). Symmetries and quantum numbers 
of BB states are explained in detail in In this work 
we focus on the two attractive channels between ground 
state static- light mesons (B mesons): a scalar isosinglet 
(cf. Fig. [Ha)), which is most attractive, and a vector 
isotriplet (cf. Fig. QJb)), which is less attractive. 

To motivate an ansatz to fit the lattice results for 
these potentials, first remember that the pair of heavy 
antiquarks is immersed in a cloud of two light quarks. 
The size of this cloud is crucial for the heavy antiquark- 
antiquark interaction. When the two antiquarks are 
much closer than twice the typical light quark cloud ra- 
dius in a heavy-light meson, the antiquark-antiquark in- 
teraction is a typical diquark interaction (cf. Fig. HJa)). 
However, when the separation of the two antiquarks be- 
comes larger, the light quarks screen the heavy antiquark 
charges as in Fig. [^b) and the antiquark-antiquark in- 
teraction vanishes. 



(a) scalar isosinglet: a = 0.29 + 0.03, p = 2.7 + 1 .2, d/a = 4.5 + 0.5 
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(b) vector isotriplet: cc = 0.20 ± 0.08, d/a = 2.5 + 0.7 (p = 2.0 fixed) 




Figure 1: (Color online). The static antiquark-antiquark po- 
tential as a function of the separation in units of the lattice 
spacing a ~ 0.079 fm. (a) Three parameter fit (a, d, p) of 
ansatz JT} to the most attractive channel, the scalar isosin- 
glet. (b) Two parameter fit (a, d; p — 2.0 fixed) of ansatz ([lj 
to the less attractive vector isotriplet. 
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Figure 2: (Color online). Screening of the antidiquark flux 
tube interaction. We show two scenarios: (a) the heavy anti- 
quarks are close and the screening of the light quarks has little 
effect on the antiquark-antiquark interaction; (b) the heavy 
quarks are well separated and their color charge is totally 
screened by the light quark wavefunctions. 



The diquark interaction for systems containing only 
static quarks has also been studied by lattice QCD. 
Flux tubes have been observed for static baryons, static 
tetraquarks and static pentaquarks [ill, E3 ■ It seems es- 
tablished that at large separations r the diquark potential 
is linear and confining ~ or. At small separations r the 
potential is rather Coulomb- like, i.e. ~ —a/r. A simi- 
lar system, which has been studied in even more detail 
by lattice methods is the "static-static meson" or static 
potential; here a sa 0.3. . .0.4, while a w (0.44 GcV) 2 is 
estimated from quark model fits and often used to set the 
scale. For two heavy antidiquarks in a cloud of two light 
quarks we expect a similar Coulomb-like potential of or- 
der —a/r at small separations. At larger separations (in 
contrast to the purely static case) the potential should be 
screened by the light quarks, as discussed in the previous 
paragraph and illustrated in Fig. [2jb). 

The screening of the heavy color charge is due to the 
decrease of the wave function ip of the light quark with 
respect to its separation from the heavy quark. One ex- 
pects this decrease to follow an exponential of a power 
law, i.e. ip(r) oc exp(— (r/d) p ), where d characterizes the 
size of the quark-antiquark system, i.e. a B meson. If the 
quark-antiquark interaction inside a B meson is domi- 
nated by a Coulomb-like term, the wave function is sim- 
ilar to that of a hydrogen atom, i.e. p = 1. In case the 
potential is rather linear, the non-relativistic Schrodingcr 
equation is solved by Airy functions corresponding to 
p = 3/2. A similar, but relativistic treatment of the 
light quark yields p = 2 instead. 

The above considerations suggest the following ansatz 
to model the heavy antiquark-antiquark potential: 



V(r) 



a 



■ exp 



r\P 
d 



(1) 



where it is expected that a«0.3...0.4, d is around half 
the size of a B meson, i.e. d < 0.5 fm and p w 1.0 . . . 2.0. 

We perform \ 2 minimizing fits of the ansatz dTJ) to 
the lattice results for the heavy antiquark-antiquark in- 
teractions shown in Fig. [T] Notice that data points for 
separations r/a = and r/a = 1 are excluded from the 
fits, because they suffer from lattice discretization errors. 
For r/a > 2 it has been checked, that these discretization 
errors are negligible compared to the statistical errors by 
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using different static quark actions |8j. On the other 
hand, data points for large separations have little rele- 
vance since the potential vanishes exponentially fast, due 
to screening, while the statistical error remains similar. 
Thus we utilize the points at distances r/a = 2, 3,4, 5, 6 
for our fits. 

For the scalar isosinglet we are able to determine all 
three parameters a, d and p via fitting. For the vector 
isotriplet a three parameter fit is not stable; therefore, 
we only fit two parameters, a and d, while fixing the ex- 
ponent to its expected value p = 2.0. The fits are also 
shown in Fig. [TJ while numerical results are collected in 
Table U (statistical errors for a, d and p have been de- 
termined via an elaborate Jackknife analysis starting on 
the level of the bare lattice correlation functions). The 
uncorrelated x 2 /dof <C 1 indicates consistency of the lat- 
tice data and our ansatz ([1]). Moreover, note that the 
resulting values for a, d and p are in agreement with 
phcnomenological expectations. 

Table I: x 2 minimizing fit results of the ansatz flj to the 
lattice static antiquark-antiquark potential; fitting range 2 < 
r/a < 6; lattice spacing a ~ 0.079 fm 



channel 


a 


d/a 


P 


X 2 /dof 


scalar isosinglet 


0.293(33) 


4.51(54) 


2.74(1.20) 


0.35 


vector isotriplet 


0.201(77) 


2.48(69) 


2.0 (fixed) 


0.06 



The potential of Eq. ([1} with the fit parameters from 
Table U corresponds to the energy of a static-static-light- 
light four quark system minus the energy of a pair of 
static-light ground state mesons. To obtain the energy 
of a heavy- heavy- light-light system, where the antiquarks 
have a heavy, but finite mass, one also needs to consider 
a kinetic term for the heavy antiquarks resulting in the 
Hamiltonian 



H = £- + 2m B + V(r), 



(2) 



where \x is the reduced antiquark mass. 

Notice that, because of screening, at large separations 
each heavy antiquark carries the mass of a B meson and 
thus [i = nifl/2, whereas at small separations it carries 
just the energy of a heavy quark [i = m^jl. When in- 
vestigating the existence of four quark bound states, we 
consider both mass values, which differ by around 5% 
(mp = 5279 MeV @], m b = 4977 MeV in quark models 
[l3[). Another possible source of systematic error is as- 
sociated with the physical value of the lattice spacing. 
This error is introduced, when converting the "size pa- 
rameter" d from dimensionlcss lattice units to physical 
units. We investigate the magnitude of this error, by us- 
ing values for the lattice spacing determined by rather 
different scale setting procedures: a = 0.079 fm is used in 
many ETMC publications and is obtained from m n , f„ 
and chiral perturbation theory [14j, while a = 0.096 fm 
corresponds to identifying the lattice result for the Som- 
mcr parameter rp with its typical value 0.5 fm. 



To get an analytical qualitative understanding, we 
first derive an approximate rule for the existence/non- 
existence of a bound state applying the Bohr-Sommerfeld 
relation corrected by the WKB approximation, to include 
the zero point energy of 1/2. For the radial equation in 
three dimensions we get, 



(3) 



where C\ and are the classical turning points. Insert- 
ing (|T|), specializing to angular momentum I = and 
integrating, yields the condition for having at least one 
bound state, 



jiad > 



9^ 2 



128 x 2 1 /pr 2 (l + l/2p)' 



(4) 



The right hand side of Eq. (j4]) has a rather moderate 
dependence on the exponent p. For example, when p in- 
creases from the expected values of 1.0 to 2.0, the right 
hand side only changes from 0.44 to 0.60. Thus the exis- 
tence of a bound state mainly depends on the product of 
parameters \xad. With the fit parameters for the scalar 
isosinglet from Table U and a bottom quark one roughly 
obtains pad ~ 2.5 GeV x 0.29 x 4.5 x 0.079 fm « 1.3 > 
0.44 . . . 0.60, which is a strong indication for the existence 
of a bound state. A similar calculation for the vector 
isotriplet yields fiad « 0.5, i.e. the existence of a bound 
state in this channel is rather questionable. 

Note, however, that these are only crude analytical 
estimates, because e.g. the WKB approximation is ques- 
tionable, when the potential is divergent at the origin 
as for a Coulomb-like potential. To investigate the exis- 
tence of a bound state rigorously, we numerically solve 
the Schrodingcr equation with the Hamiltonian The 
strongest binding is expected in an s-wave, for which the 
radial equation is 



' 2[i dr 2 



+ 2m B + V(r) 



R{r) = ER(r) (5) 



with the wave function ip = ip(r) = R(r)/r. We im- 
pose Dirichlet boundary conditions i?(r max ) = at suffi- 
ciently large r max (we checked that results are stable for 
r max > 10 fm). The radial equation §5§ can be solved by 
standard methods (e.g. 4th order Runge-Kutta shooting) 
up to arbitrary numerical precision. 

In accordance with our analytical estimates we find 
binding for heavy bottom-bottom tetraquarks in the 
scalar isosinglet channel. The binding energy E — 2ms ~ 
30 MeV ... 57 MeV depends to some extent on the re- 
duced mass (either /j = m^/2 or fi = ttib/2) and on 
the value of the lattice spacing (either a = 0.079 fm or 
a = 0.096 fm); cf. Table [TTl Note that the binding energy 
is negative by around 1.76cr . . . 3.00c This confirms the 
existence of an heavy-heavy-light-light tetraquark with a 
confidence level of >96%. The probability to find the 
two heavy antiquarks at separation r is proportional to 



4 



Table II: Binding energy and confidence for the existence of 
a heavy-heavy-light-light tetraquark for the scalar isosinglet 
channel. 



M 


a 


E - 2m B 


confidence for binding 


m b /2 


0.079 fm 


-30(17) MeV 


1.76 cr, 96% binding 


m b /2 


0.096 fm 


-49(17) MeV 


2.88 cr, 100% binding 


ma/2 


0.079 fm 


-38(18) MeV 


2.11 cr, 98% binding 


ros/2 


0.096 fm 


-57(19) MeV 


3.00 cr, 100% binding 



probability to find the b antiquark pair at separation r 





fi = m b /2, a = 0.079 fm 




\i = mtj/2, a = 0.096 fm - 




H = m B /2. a = 0.079 fm 




s£\ n = m B /2. a = 0.096 fm 







0.6 
r in fm 



Figure 3: (Color online) . Radial probability density oc \R(r) 2 
for the separation of the heavy antiquarks (scalar isosinglet). 



Further possible systematic errors are not expected 
to weaken the binding. The static approximation of a 
ground state with bottom antiquarks is valid, since the 
bottom quark mass is a very hard scale compared with 
the scale Ang p. Moreover, lattice computations of B 
mesons [HI, Hg an d b baryons [l7| within the same lat- 
tice setup showed that static quarks are a rather good 
approximation of bottom quarks. Finite volume effects 
are exponentially suppressed in lattice results and the 
spatial extension of the lattice L w 1.9 fm seems large 
compared to the typical size of the tetraquark, which is 
related to d m 4.51 a rj 0.36 fm and the average separa- 
tion of the heavy antiquarks, which is w 0.25 fm. Finally, 
the light u/d quark masses in the lattice computation are 
unphysically heavy (corresponding to w 340 MeV) . 
Decreasing the light quark masses to their physical val- 
ues should increase the light cloud radius of a heavy-light 
meson and, therefore, lead to stronger binding (cf. Fig. 
Ha)). 

To summarize, we find clearest evidence for the exis- 
tence of an antibottom-antibottom-light-light tetraquark 
bound state. This result is very promising with respect 
to further lattice QCD calculations with higher statistics 
or even lighter dynamical quarks, and for experimental 
searches of doubly bottom hadrons. 



|i?(r)| 2 and shown in Fig. [31 The average separation is 
around 0.25 fm. 

Solving Schrodinger's equation (|5|) for the vector 
isotriplet gives strong indication that the potential of 
this channel is too weak to generate a bound state, i.e. 
a tetraquark. The resulting wave function is essentially 
a plane wave with positive energy E — 2m b within more 
than 10 cr. 
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